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Abstract 

Lattice QCD calculation of s-wave hadron scattering lengths in the channels 
TT-TT, T^-N, K-N, K-N and N-N is carried out in the quenched approximation at 
i5 — ^1 — 5.7. A variant of the method of wall source is developed for this 
purpose, which reduces the computer time by a factor on an x T lattice 
compared to the conventional point source method and avoids the Fierz mixing 
problem. A version of the method in which gauge configurations are not fixed to 
any gauge can be extended to calculate disconnected quark loop contributions in 
hadron two- and three-point functions. An analytical estimate of statistical errors 
for this method is worked out, and the magnitude of errors without and with 
gauge fixing is compared for the case of tt-tt four-point functions calculated with 
the Kogut-Susskind quark action. For tt-tt scattering both 7 = and 2 scattering 
lengths are evaluated using the Kogut-Susskind and Wilson quark actions on a 
12^ X 20 lattice. For the same size of lattice, tt-A^, K-N and K-N scattering lenghts 
are calculated with the Wilson quark action. For the tt-tt and tt-A^ cases simulation 
results are consistent with the predictions of current algebra and PCAC within 
one to two standard deviations up to quite heavy quark masses corresponding 
to rriT^/mp 0.74, while for the K-N and K-N cases the agreement is within 
a factor of two. For A^-A^ scattering a phenomenological study with one-boson 
exchange potentials indicate that the deuteron becomes unbound if the quark mass 
is increased beyond 30-40% of the physical value. Simulations with the Wilson 
action on a 20^ lattice with heavy quarks with m^/mp ^ 0.74 — 0.95 show that the 
nucleon-nucleon force is attractive for both spin triplet and singlet channels, and 
that the scattering lengths are substantially larger compared to those for the tt-tt 
and TT-A" cases even for such heavy quarks. Problem of statistical errors which has 
to be overcome toward a more realistic calculation of hadron scattering lengths is 
discussed. 



2 



1 Introduction 



Calculation of scattering lengths of hadrons presents a variety of problems of 
physical interest in lattice QCD approach to strong interactions. For the case of 
TT-TT and TT-N scattering, the experimental s-wave scattering lengths are of order 
0.1 — 0.4fm which are quite small compared to the sizes of the pion and nucleon 
(see Table 1 for experimental values). It is well known that the small values can 
be understood as a result of soft pion theorems [2] that follow from approximate 
chiral symmetry of up and down quarks realized in a spontaneously broken mode; 
scattering amplitudes involving pions vanish at threshold for zero pion mass, and 
the s-wave scattering lengths predicted for the experimental pion mass are in 
reasonable agreement with experiment as recapitulated in the last column of Ta- 
ble 1. Whether a dynamical calculation directly based on the QCD Lagrangian 
successfully explains experiment therefore provides a useful testing ground of our 
theoretical understanding of chiral symmetry based on lattice QCD. 

The situation is quite different for nucleon-nucleon scattering. Experimentally 
the scattering lengths are very large, being of the order of lOfm (see Table 1). 
Since chiral symmetry places no constraint on the low energy behavior, the large 
scattering lengths are purely dynamical phenomena. So far theoretical approach 
to the low energy nucleon-nucleon scattering has remained at a phenomenological 
level in which meson exchange models with tunable coupling constants are utilized 
to reproduce the low energy phase shifts[3]. Deriving the large scattering lengths 
from first principles of QCD represents an important challenge to lattice QCD. 

Calculation of scattering lengths poses several problems from the technical 
point of view. One of the problems is that Euclidean hadron Green's functions 
amenable to Monte Carlo evaluation methods of lattice QCD are generally only 
indirectly connected to real time scattering amplitudes. An elegant solution to 
this problem is the formula[4, 5] relating scattering phase shifts, and hence also 
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scattering lengths, to the energy of two hadron states on a finite spatial lattice 
which can be extracted from the exponential decay of Euclidean hadron four-point 
functions in time [6]. 

Another problem, which is computationally quite troublesome, is the calcula- 
tion of hadron four-point functions itself. For hadrons in a definite momentum 
state, such a calculation generally requires quark propagators connecting two ar- 
bitrary space-time sites. With the conventional method of point source the nec- 
essary number of quark matrix inversions equals the space-time lattice volume, 
which would require a prohibitively large amount of computer time. This is a 
novel situation quite different from a calculation of connected three-point func- 
tions of hadrons, which can be reduced to that of connected two-point functions 
through source methods [7]. We note that a similar difficulty arises even for two- 
or three-point functions when one attempts to calculate contributions involving 
disconnected quark loops. Typical examples are the two-quark loop amphtude for 
the fiavor singlet rl propagator and fiavor singlet nucleon matrix elements such as 
the 'K-N a term. 

We found that a variant of the wall source technique [8, 9] allows to overcome 
the problem with a modest cost of computing power[10, 11, 12, 13, 14, 15, 16]. 
Particularly interesting is a modified version of the original proposal[8] which does 
not fix gauge configurations to any gauge. This version can be applied to calculate 
disconnected quark loop contributions in two- and three-point functions as well as 
hadron four-point functions. Our chief findings obtained with the method for the 
TT-TT, TT-iV and A^-iV scattering lengths have been briefly reported in Refs. [10, 11, 
12]. In this article we present the full details of the method and analyses of our 
calculation of hadron scattering lengths. We also report additional results for K- 
N and K-N scattering. For applications of the method to calculate disconnected 
quark loop contributions in hadron two- and three-point functions, we refer to our 
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work on the r( meson mass[13], the 'k-N a term[14, 15] and the axial vector matrix 
elements of the proton[15, 16]. 

Historically the first attempt toward a lattice calculation of hadron scattering 
lengths was made by Guagnelli, Marinari and Parisi[17] for the tt-tt and 'k-N cases 
using the Kogut-Susskind and Wilson quark actions in quenched QCD. They found 
a finite volume shift of the energy of two-hadron states and examined the volume 
dependence. However, they calculated only gluon exchange diagrams, and thus 
the results cannot be compared with physical scattering lengths. Gupta, Kilcup, 
Patel and Sharpe[18, 19] developed an analysis of the physical tt-tt scattering length 
for the 1 — 2 channel which can be calculated with the conventional wall source 
method. Their studies were made in quenched QCD using the Kogut-Susskind[18] 
and Wilson [19] quark actions, and lattice results for scattering lengths were found 
to be consistent with the predictions of current algebra and PCAC for the 1 — 2 
channel up to quite large quark masses. 

The focus of the present work consists in an extension of the previous 1 — 2 
results for the tt-tt scattering length to the 7 = channel and also to the t^-N 
(/ = 1/2,3/2), K-N and ^-A^ (/ = 0, 1) scattering. These extensions require the 
modified version of the wall source method developed here. We also carry out an 
exploratory study of the nucleon-nucleon scattering lengths. Our simulations have 
been made within quenched QCD at /9 = 5.7 on a 12^ X 20 lattice for the TT-TT, 'K-N , 
K-N and K-N cases, and on a 20^ lattice for the N-N case. For tt-tt scattering 
we employed both the Kogut-Susskind and the Wilson quark action, while other 
scattering lengths are studied with the Wilson quark action. 

This paper is organized as follows. In Sec. 2 we describe the formalism for 
calculation of hadron scattering lengths including the Liischer's formula [5] and 
our calculational technique of modified wall sources for construction of hadron 
four-point functions. In Sec. 3 we summarize our data sets. The results for tt-tt 
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scattering lengths are presented in Sec. 4, where we also compare Kogut-Susskind 
results with and without gauge fixing, and discuss the question of infrared sin- 
gularities in quenched QCD associated with the degeneracy of 77' and tt. Results 
for TT-A^, K-N and K-N scattering lengths are given in Sec. 5, while those for 
N-N scattering lengths are presented in Sec. 6. In Sec. 7 we discuss the issue of 
statistical errors in relation to a fully realistic lattice QCD calculation of hadron 
scattering lengths in the future. Our conclusions are summarized in Sec. 8. 

Throughout the article we employ lattice units for expressing physical quanti- 
ties and suppress the lattice spacing a unless necessary. 

2 Formalism and method of measurement 

2.1 Scattering length from two-particle energy in a finite 
box 

The basic formula for a calculation of scattering lengths was derived by Liischer[5] 
who related the s-wave scattering length ao between the two hadrons hi and /i2 to 
the energy shift of the two hadron state at zero relative momentum confined in a 
finite periodic spatial box of a size . The formula is given by 

Ehih2 - [rrihi + mh^) ^ 7^ (1 + Ci— + C2(— ) ) + 0(L ) (1) 

with ci = -2.837297, C2 = 6.375183. 

The energy E^^h^ is extracted from the large time behavior of the hadron four- 
point function with the hadrons in the zero momentum state defined by 

X4 X3 X2 S\ 

(2) 

In order to avoid Fierz rearrangement of quark lines to be discussed in Sec. 2.4, 
we choose — t -\-l,t2, — t,t2 — 1 and ti — 0. In this case the large time behavior 

reads 

Ch^h2{t + 1, 1, 0) = Zhih2^xp{-Eh^h2t) H , (3) 
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where dotted terms denote contributions of higher excited states. 

The hadron masses and m/jj needed to obtain the energy shift 6E = Eh^h^ — 
{rrihi + ruh^) are calculated in the usual manner from the two-point functions, 

Cft,(t,0) =< ^O;,,(f3,t)E^Ii(^i'0) >= Zh,eM-mh,t) + ■ ■ ■ , (4) 

XZ XI 

Ch,{t + 1, 1) =< 5: Oh, (f4, ^ + 1) E < (^2, 1) >= Zh,eM-mh,t) + • • • . (5) 

X4, X2 

We extract the energy shift SE from the ratio: 

""^'^^ c,,mc,,(t + 1,1) ^z^r"^-'^'^^-"- 

2.2 Hadron four-point functions 

Consider s-wave tt-tt scattering, for which isospin and 2 channels are allowed 
due to Bose symmetry. We construct the two-pion operators for these isospin 
eigenchannels as 

oZoit) = ^ {o^At)o.-{t + 1) - oAt)OAt + 1) + 0^-{t)0^+{t + 1)) , (7) 

OZ2{t) = a+Wa+(i + l), (8) 
with the pion operators defined by 

On+{t) = -Y.d{x,th5u{x,t), (9) 

X 

0^-{t) = ^u{x,t)^,d{x,t), (10) 

X 

O^o{t) = (^(^' ^h^uix, t) - d{x, t)-f5d{x, t)) . (11) 

In Fig. 1 we display quark line diagrams contributing to the tt-tt four-point func- 
tions, denoting them as direct (D), crossed (C), rectangular (R) and vacuum (V) 
diagrams (meaning of circles and bars will be explained in Sec. 2.3). The ratio 
R{t) for isospin eigenchannels can be expressed in terms of the four amplitudes in 
the following combination; 

RZoit) = i?^(t) + V(t)-3i?«(t) + ^i?^(t), (12) 
RZ2{t) = i?^(i)-i?^(i). (13) 



For T^-N scattering the 7=1/2 and 3/2 operators may be taken as 

Oti/2{t) = ^OMt+l)Or.{t)--^OAt+mp{t), (14) 
0]^,/,{t) = 0^+{t + l)0,{t). (15) 

For the pion operators we use (9-11). Relativistic operators for the proton (p) 
and neutron (n) are defined with the charge conjugation matrix C = 7472 taking 
account of the symmetries for isospin and non-relativistic spin representations, 

Op{t) = J2^abc [{'u%X,t)C-'^^d\x,t)) U\x,t) - ('d%X,t)C-^^5u\x,t)) u'{x,t) 

X 

(16) 

Onit) = Y^Sabc [(V(f,i)C-S5f^'(^,i)) d''ix,t) - {^d''{x,t)C-^^^u\x,t)) d%X,t) 

X 

(17) 

In terms of O'^^ we construct the t:-N four-point function according to 



(18) 



where trace refers to Dirac components. 

Topologically four types of diagrams contribute to n-N four-point functions, 
direct (D), crossed (C), rectangular (R) and crossed rectangular (CR), as shown 
in Fig. 2. In contrast to the tt-tt case, these diagrams have 6, 18, 18 and 36 
members having the same quark line topology but different quark contractions. 
All four types of diagrams contribute to both isospin eigenchannels. Weights of 
the members are listed in Appendix. 

For 7 = and 1 operators for K-N and K-N scattering we take 

OfZit) = ^{OK+{t+l)On{t)-OK'^{t + l)Op{t)), (19) 

Ofi(t) = OK+{t + l)0,{t), (20) 

ofm = -^{0^o{t + l)On{t)-OK-{t + ^)OM. (21) 
O0(t) = O^o(t + l)0p(t), (22) 



where the nucleon operators Op and C„ are given in (16-17) and the K meson 
operators are defined by 

OK+{t) = E^(^'^)75«(^,^), (23) 
OMt) = X^s(f,i)75d(f,i), (24) 
OKo{t) = -^ci(f,i)75s(f,i), (25) 

X 

OK-{t) = Y.^{x,t)^,s{x,t). (26) 

X 

In contrast to the n-N case, K-N scattering receives contributions only from 
direct (D) and crossed (C) type diagrams and K-N scattering only from direct 
and rectangular (R) ones. Weights of various quark contractions to isospin eigen 
amplitudes are given in Appendix. 

In s-wave N-N scattering spin triplet ^5*1 and singlet ^5*0 states choose isospin 
and 1 respectively due to Fermi statistics. For these eigenchannels we use the 
operators, 

O^S^m = l=^Op{t)On{t + l)-On{t)Op{t + l)), (27) 
OCSo){t) = Op{t)Op{t + l), (28) 

where the upper two Dirac components of the nucleon operators Op and On as 
defined in (16-17) are combined to form the respective spin states. N-N four- 
point functions require only two types of diagrams shown in Fig. 3. The direct (D) 
diagram has 36 members with different quark contractions. For the crossed (C) di- 
agram, which has 324 members, diagrams with a single quark-antiquark exchange 
and a double quark-antiquark exchange are equivalent because the two nucleon 
operators are summed over all spatial sites for s-wave states. Since quark contrac- 
tions for N-N four-point functions are very complicated and tedious, we employed 
symbolic manipulations on computers to work out the weights of amplitudes. 
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2.3 Wall source method without gauge fixing 

Let us consider the numerical procedure for a calculation of hadron four-point 
functions. For tt-tt scattering shown in Fig. 1 the direct and crossed diagrams 
can be easily calculated because we need only two quark propagators with wall 
sources placed at the fixed time shces ti and i2[18, 19] in order to construct the 
corresponding four-point amplitudes for arbitrary values of ts and (the A'"-A'^ case 
is similar) . This does not apply to the case of rectangular and vacuum diagrams 
which require additional quark propagators connecting the time slices ^3 and ^4. 
The same difficulty also exists in the n-N, K-N and K-N cases for the rectangular 
and crossed-rectangular types of diagrams. 

We handle this problem by calculating T quark propagators on an x T 
lattice, each propagator corresponding to a wall source placed at the time slice 
t = 0, • • • , T — 1, which are defined by 

Y,Dn',n"Gt{n") = E5n',(a-,t), < t < T - 1 (29) 

n" X 

where D denotes the quark matrix for the Wilson or the Kogut-Susskind quark 
action. The combination of Gt{n) that we employ for hadron four-point functions 
are displayed in Figs. 1-3, where short bars represent the position of wall source and 
circles the sink. For example the tt-tt rectangular diagram in Fig. 1(c) corresponds 
to 

C (t4, ts, ti) = 

5^ < ReTr[Gl,(f2,t2)G't,(f2,t2)G'Uf3,t3)Gti(5^3,t3)] >, ^^^^ 

where daggers mean conjugation by 75 for the Wilson quark action and by the 
even-odd parity (—1)" for the Kogut-Susskind quark action. 

Using the relation Gt{n") — J2xD~,} (^(), we see that this prescription yields 
four-point amplitudes corresponding to non-local and non-gauge invariant hadron 
operators at the time slices where two wall sources or one wall and one sink are 
placed {e.g., at the time slices ti and in (30)). These terms create gauge- variant 
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noise. One way to suppress the noise is to fix gauge configurations to some gauge 
as is done in all recent work using wall sources [9]. A potential drawback is that 
gauge non-invariant states may contaminate the four-point function. 

Alternatively one can choose not to fix gauge configurations to any gauge since 
gauge dependent fluctuations should cancel out in the ensemble average. This 
is the wall source version of the original proposal of extended sources [8]. One 
might worry if signal stands out among noise. For example, if two wall sources are 
placed at the same time slice, there are 0((L^)^) gauge dependent non-local terms 
relative to 0{L^) local gauge invariant ones, and hence naively the magnitude 
of the noise 0{^J{L'^Y) would be comparable to that of the signal. However, a 
generalization[20] of the well-known argument[21] can be used to show that the 
noise is in fact smaller than signal by a factor L^/^ for sufficiently large L. 

Let us illustrate this point with a simpler example of a two-point function Cr(t) 
for the meson operator qVq for the Wilson quark action. We evaluate this quantity 
as 

Cvit) = < E Tr {gI{x, t)TG^{x, t)T) >, (31) 

where the factor 1/L^ is inserted to ensure a finite limit as L — > oo and the trace 
is taken over Dirac indeces and color indeces. Substituting the expression 

G,{x,t)^Y.G{{x,t),{y,Q)) (32) 
y 

for the wall-source quark propagator GQ{n) in terms of the point-to-point quark 

propagator G{n,m), we can rewrite (31) as 

<^r(t)=<ETr(Gr(^,^)r)>, (33) 

where 

Gr{y, z)^^Y.G {{y, 0), (x, t)) TG {{x, t), {z, 0)) . (34) 

X 
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The dispersion of the two-point function Cr{t) is given by 

a' = <|ETr(Gr(^,i')r)r>-|<ETr(Gr(f,i')r)>r 
y,^ y,z 

= < E ^ (Gr(y, m E Tr {Gr{z', y')r) > 

y,z z',y' 

- < J2^HGr{y,m X E Tr(Gr(i",y')r) > (35) 

y,z z',y' 

where we have used the relation G'^(n, m) = 75G(m, 71)75 and assumed that 
75r75 = ±r. We recognize the first term on the right-hand side to be the ex- 
pression for the four-point function with a quark operator placed at z and y' and 
an anti-quark operator at y and z' of the time slice t = 0. Since gauge configura- 
tions are not fixed to any gauge, non- vanishing results are obtained only ii z — y 
and y' — z' or z — z' and y' = y. We thus find that 

(7^ = 4 + 4, (36) 

with 

al = ^<Tr(Gr(f,^)r)Tr(Gr(i',i')r)> 

y,z 

- ^ < IV {Gr{y, m X Tr {Gr{z, z)r) > (37) 

y,z 

al = 5:<Tr(Gr(y,i')r)Tr(Gr(^^,y)r)> 
y,z 

= 1^ E E<Tr(Gr(y,^)mGr(^-',y-')m)> (38) 

-■-^ i=S,V,T,A,P y,z 

where for we used the Fierz transformation for Dirac and color indeces to com- 
bine the two traces into a single trace and ignored color non-singlet combinations. 
Substituting (34) it is easy to see that ajj equals the direct amphtude for the 
meson-meson scattering in the channel F -|- F — > F -|- F with all meson at rest, mul- 
tiplied by a factor Similarly equals the crossed amplitude for the channel 
rVi + rVi -^r + r multlpfied by We thus find that the statistical error for 
the two-point function evaluated with Nconf independent gauge configurations is 
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given by 

SCr{t) 



conf 



\ AT 7-3 I ^r+r-.r+rW + t;^ X] ^rri+rri-^r+rW ) ■ (39) 

\^ ^^conf-l^ y J-^ i=S,V,T,A,P ) 

This formula immediately shows that the statistical error decreases as L~^/^ 
relative to the signal. The direct amplitude decreases as exp(— 2mri) up to poly- 
nomials in t with mr the mass of the meson in the channel F. For the crossed 
amplitude we expect the dominant term in the sum over % to come from the chan- 
nel TT -I- TT — > F -I- F having the lightest energy. For a sufficiently large spatial size L 
the crossed amplitude in this channel may be approximately calculated by a con- 
volution of free propagation of two pions from the time slice i = to i'(0 < < i), 
a local ttttFF coupling dXt — t' and free propagation of two mesons in the channel 
F from t = t' to t. As an estimate of the relative error we then find 

It is straightforward to generalize the argument to hadron four-point functions 
Chihi- The relative error is given by 



5Ch-j,,{t) ^ 1 exp((m,,, + TUh, - {uh, + nh,)m^)t). (41) 

where Uh = 1 for mesons and 3/2 for baryons. 

We should note that the magnitude of the proportionality constant in (40) and 
(41) is not known, and hence the magnitude of errors may vary depending on the 
channel. With a practically manageable statistics of a few hundred configurations, 
we found that statistical errors are small for pion two- and four-point functions, 
whereas signal deteriorates for p and nucleon. For nucleon, in particular, errors are 
significant even for the propagator. We, therefore, use the Coulomb gauge fixing 
for the nucleon source for calculation of scattering lengths involving nucleon. 
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2.4 Avoidance of Fierz contaminations 

Our method of calculating quark propagators for all possible temporal posi- 
tions of wall sources has another advantage. Suppose that we ony have the quark 
propagator corresponding to a single wall source placed at the time slice t. Cal- 
culation of meson four-point functions in this case necessarily involve two meson 
operators Z]x,y,a ^(^) t)" and Z)x',y',fe ^)*r'g(y', t)^ placed at the same 

time slice t. It is easy to see that there are two contributions, one in which each 
of the pairs {q{x,tY,q{y,tY) and {q{x' ,1)^ ,q{y' ,1)^) form color singlets, and the 
other in which the Fierz-rearranged pairs t)", q{y\ t^) and {q{x', tY, q{y, t)") 
form color singlets. Such a rearrangement causes a mixing of diagrams having dif- 
ferent quark line topologies. In general it also switches the spin-parity of the two 
hadrons in the initial state. This leads to quite a complication of analyses, espe- 
cially for the Kogut-Susskind quark action, as was discussed in detail in Ref. [18]. 
Disentangling the two contributions is possible in principle, but difficult in prac- 
tice. This, in fact, is the reason why the work of Ref. [18] could not determine the 
direct and crossed amplitudes separately for the tt-tt four-point function for the 
Kogut-Susskind quark action. 

The problem can be trivially solved in our method since we calculate quark 
propagators for all possible temporal positions of wall source, and therefore hadron 
operators can be placed at different time slices. In practice we separate hadron 
operators by a unit time slice both for the initial and final pairs of hadrons. 

3 Data sets 

In Tabic 2 we list the parameters of our simulation. All of our calculations 
are made in the quenched approximation at /9 = 5.7. The inverse lattice spacing 
determined from the p meson mass in the chiral limit equals — 1.44(2)GeV 
for the Wilson quark action and = 0.98(ll)GeV for the Kogut-Susskind quark 
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action. We mostly employed a 12^ x 20 lattice for tt-tt, tt-A^, K-N and K-N scat- 
tering and a 20"^ lattice for N-N scattering, anticipating large scattering lengths 
for the latter. The lattice size must be large enough so that we can employ as weak 
couphng as possible to avoid finite lattice spacing effects, yet it should not be too 
large so as not to spoil a detection of a small energy shift of 0{L~^) predicted 
by the Liischer's relation (1); our parameters are a reasonable compromise with 
the present computing resources. Gauge configurations are generated with a 5-hit 
pseudo heat bath algorithm, discarding 1000 configurations for thermalization and 
employing every 1000th configuration for hadron four-point function analyses. 

For the tt-tt four-point function we made bulk of our calculations employing wall 
sources without gauge fixing. A calculation with Coulomb gauge fixing was also 
carried out for a subset of configurations with the Kogut-Susskind quark action in 
order to compare the two methods (see Sec. 4.1.2). In the case of t:-N scattering 
the pion source is treated with the wall source method without gauge fixing, while 
for the nucleon source placed at i = we fixed the t — Q time slice of gauge 
configurations to Coulomb gauge in order to enhance signal to noise ratio. The 
K-N and ^-A^ four-point functions are calculated in the same way. We used the 
same quark mass for strange and up-down quarks for this case. For calculations 
of the N-N four-point functions gauge configurations are fixed to Coulomb gauge 
over the entire space-time lattice as this leads to better signals. 

Quark propagators are calculated with the Dirichlet boundary condition in 
the temporal direction in order to avoid contamination from hadrons propagat- 
ing backward in time. The periodic boundary condition is used in the spatial 
directions. For the Kogut-Susskind quark action the standard conjugate gradient 
algorithm is employed for inverting the quark matrix, and for Wilson quark ac- 
tion the ILUCR method[22]. The stopping condition is chosen to ensure a 0.1% 
accuracy in hadron four-point functions for each configuration. 
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In Table 3 we summarize tt, p and nucleon (iV) masses in lattice units ob- 
tained on our ensemble of configurations through standard single exponential fits 
of hadron propagators over the time interval 6<f<12 — 14. Also listed are the 
values for the pion decay constant obtained with the tadpole-improved renormal- 
ization factor[25] using otjig{l/a) = 0.2207 for the coupling constant. The values 
listed are in agreement with the results obtained on larger lattices[23, 24], which 
are also given in Table 3. A necessary condition for the apphcability of the formula 
(1) is that the lattice size L is sufficiently large so that finite-size effects for hadron 
masses are negligible. The agreement provides a check, albeit indirect, on this 
point. 

We estimate errors by the single elimination jack-knife procedure for all quan- 
tities including hadron masses, energy shifts and scattering lengths obtained by 
fits of two- and four-point functions. 

4 7V-7T scattering lengths 

4.1 Kogut-Susskind quark action 

4.1.1 No gauge fixing 

Current algebra and PCAC predict [2] that the s-wave tt-tt scattering lengths, 
to leading order in m^, take the values. 



7=0 


7 rriTT 


'0 — 


^327r /2 


7=2 _ 


2 


'0 


327r /2 



It has been shown in Ref. [18] that these results can also be derived for the Kogut- 
Susskind quark action at a finite lattice spacing using Ward identities for U{1) chi- 
ral symmetry under some continuity assumptions on the pion four-point function 
in external momenta. Whether numerical simulations yield results in agreement 
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with (42-43) therefore provides a valuable check of the lattice method for hadron 
four-point functions. 

The use of Kogut-Susskind action has a subtle problem concerning the inter- 
pretation of the four quark flavors corresponding to each Kogut-Susskind field. 
Following Ref. [18] we introduce a Kogut-Susskind field for each continuum flavor, 
regarding the four flavors of the Kogut-Susskind action as a four-fold duplication 
of each continuum flavor. In this view the expression (12-13) for the tt-tt four-point 
function in isospin eigenchannels needs to be modified to 



where the factors Nf = A compensate a different number of Kogut-Susskind fiavor 
traces in the four types of diagrams. For the pion operator it is most natural 
to take the one in the Nambu-Goldstone channel corresponding to C/(l) chiral 
symmetry. This is the choice for which the current algebra result (42-43) can be 
derived. 

We have carried out simulations on an 8^ X 20 lattice and a 12=^ x 20 lattice, both 
at /? = 5.7 and niq = 0.01 in quenched QCD, employing the method of wall source 
without gauge fixing. The results on a 12^ x 20 lattice have been briefiy reported in 
Refs. [10, 11]. In Fig. 4 we recapitulate the individual ratios R^{t) {X — D,C,R 
and V) on a 12'^ x 20 lattice plotted as functions of t. Good signals observed up 
to i ~ 12 for the rectangular amplitude and up to i ~ 8 for the vacuum amplitude 
demonstrate the practical applicability of the method of wall source without gauge 
fixing. For the direct and crossed amplitudes previous results[18] employed a single 
wall source at t — 0, and could not resolve the Fierz-rearranged mixing between 
the two amplitudes. As we discussed in Sec. 2.4 separate calculation of these 
diagrams are possible in our case since the two pion sources are placed at a unit 
lattice spacing apart in the time direction. 



AT O 

Rlloit) = R''{t) + ^R^{t)-^NfR^{t) + -B^{t), 
RT=2{t) = R''{t)-NfR^{t). 



(44) 



(45) 
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Physically the interesting features in Fig. 4 are (1) a very fiat behavior of 
the direct amplitude showing that the interaction is weak in this channel, (2) 
an almost linear t dependence of the crossed and rectangular amplitudes with a 
slope of similar magnitude but opposite sign, and (3) a small value of the vacuum 
amplitude. These characteristics are in agreement with expectations from chiral 
perturbation theory and the empirical OZI rule. We also note that the crossed 
and rectangular amplitudes have the same value at t = as expected from the 
fact that their analytic expressions are identical at this value of t. 

In Fig. 5 we plot the ratio Ri{t) projected onto the 7 = and 2 isospin channels 
for a 12^ X 20 lattice. A decrease of the ratio for the 1 — 2 channel corresponds to 
a positive energy shift and hence to a repulsive interaction in this channel, while 
an increase of Ri=o{t) means attraction for the 7 = channel. A dip observed at 
t ^ 2 for the 7 = channel becomes more pronounced on an 8^ x 20 lattice as 
shown in Fig. 6(note the change of vertical scale from Fig. 5). The physical origin 
of the dip and its size dependence is not clear to us. 

Extraction of the energy shift SEj from Ri{t) and thence the s-wave scatter- 
ing length require some care. Since our calculations are made in quenched QCD, 
rescattering effects that require sea quark loops are not properly taken into ac- 
count. This effect starts with terms of 0{t'^) in Ri{t). For the Kogut-Susskind 
quark action there are further complications arising from the non-dcgencracy of pi- 
ous in the Nambu-Goldstone and other channels at a finite lattice spacing. Briefly 
stating, the contribution of non-Nambu-Goldstone pions in the intermediate states 
is exponentially suppressed for large times due to their heavier masses compared to 
that of the Nambu-Golstone pion, instead of yielding terms growing with powers 
of t for the degenerate case. This affects both the relation between the ratio Ri{t) 
and the energy shift 6Ei and that between 6Ei and scattering lengths (1). A de- 
tailed analysis[18] shows that the 0{L~^) terms in the relation (1) are invalidated 
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by the effect, and that O(t^) terms in Ri{t) are not correct even in full QCD. Thus, 
a proper extraction of scattering lengths for the present case requires the spatial 
lattice size to be large enough so that the 0{L~^) terms are small; SEj should also 
be small so that there is a range of t over which Ri{t) exhibits a linear behavior. 

Based on the considerations above we fit Ri{t) with a linear form Zi{l — dEjt) 
with the fitting range chosen to be 4 < t < 9, ignoring higher order terms. The 
fitted values of SEj and the results for the s-wave scattering length ao in lattice 
units obtained using (1) are summarized in Table 4. Here we used pion mass given 
in Table 3. The errors quoted for oq are statistical only. For a 12^ x 20 lattice the 
0{L~^) term contributes 10% in the 7 = channel, although it is negligibly small 
(< 1%) for 1 = 2. The O(t^) terms in Ri{t), neglected in the above procedure, 
contribute a few percent in both isospin eigenchannels. These uncertainties are 
within the statistical errors of 16% for 7 = and 9% for 7 = 2. 

On an 8'^ x 20 lattice, the Z factor for the 7 = channel deviate severely 
from unity due to a dip of Ri=o{t) for small t, rendering the value of the extracted 
scattering length questionable. For the 7 = 2 channel, the value for ao is consistent 
with that from a 12^ x 20 lattice. 

In Table 4 we also list in brackets the values predicted by current algebra 
and PCAC (42-43) substituting the value of and the pion decay constant 
in Table 3. We observe that the simulation results are consistent with (42-43) 
within 1 — 1.5 standard deviations, which we find quite reasonable in view of the 
systematic uncertainties discussed above. Results for ao(7 = 2) in agreement with 
(43) was previously obtained [18] with the assumption that the contribution of the 
direct diagram is negligible. 

4.1.2 Coulomb gauge fixing 

We have repeated the calculation of the tt-tt four-point function, fixing gauge 
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configurations to the Coulomb gauge for all time slices for a 12^ x 20 lattice. In this 
case the wall source and sink for pion propagators in the denominator of R^it) 
have to be chosen in the same combination as for the four-point amplitude in the 
numerator for each type of diagrams X — D,C,R and V. Without this matching 
the Z factors for the ratio R^it) would be different between the diagrams, and 
the combination (44-45) would no longer project out isospin eigenchannels. This 
problem is absent for wall sources without gauge fixing because non-local terms in 
the wall source operator cancel out after averaging over gauge configurations. 

In Fig. 7 we compare the ratio Ri{t) calculated in Coulomb gauge (open sym- 
bols) and without gauge fixing (filled symbols) for the same set of 60 configurations 
on a 12^ X 20 lattice. We suspect that a wiggle seen up to t ~ 4 for the / = 
channel in the Coulomb gauge calculation is due to a contribution of non-local p 
that can be emitted from the gauge fixed wall source. A slight discrepancy in the 
intercept for the 1 = 2 channel is ascribed to the difference of Z factors between 
the two calculations. Except for these points the two data sets are consistent. 
Results for scattering lengths extracted with the same procedure as for the case 
without gauge fixing agree with those for the case of no gauge fixing within one 
standard deviation (see Table 4). 

An interesting question is to what extent gauge non-invariant errors are con- 
trolled in the method of wall source without gauge fixing compared to the gauge 
fixed case. To examine this point we plot in Fig. 8 the single-elimination jackknife 
errors of Ri{t) for the two calculation on the same set of 60 configurations on a 
12^ X 20 lattice. Although the magnitude of errors is larger for the non-gauge fixed 
case (filled symbols) , the amount of increase of errors is contained at the level of a 
factor of 1.5 — 2 times those for the Coulomb gauge fixing (open symbols), showing 
that gauge variant noise does not give rise to a serious problem for calculation of 
TT-TT four-point functions. 
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4.2 Wilson quark action 

We have also applied the wall source method without gauge fixing to the Wilson 
quark action at the hopping parameter K — 0.164, employing 70 configurations 
on a 12^ X 20 lattice. 

In Fig. 9 we show the ratios R'^{t) for X — D,C,R and V. The direct and 
crossed amphtudes were previously calculated in Ref. [19] through gauge fixed 
wall sources placed at the same time slice. We observe in Fig. 9 that the slope for 
the rectangular amphtude R^{t) is quite small compared to that for the crossed 
amphtude R'^{t). This is different from the Kogut-Susskind case for which both 
amplitudes exhibit a similar slope (see Fig. 4). We consider that this is due to 
the heavy quark mass {mT^/mp — 0.74) for the present calculation with the Wilson 
quark action compared to the small value {rnT,/mp = 0.33) for the Kogut-Susskind 
case. Another difference between the Wilson and Kogut-Susskind results is a 
positive curvature for the direct amplitude seen in Fig. 9 for the Wilson case, 
whose origin is not clear to us. Also the errors for the vacuum diagram blows up 
much more rapidly (compare Fig. 9 with Fig. 4). This, however, may be ascribed 
to a heavier pion mass for the Wilson simulation as will be discussed in Sec. 7 
below. 

The ratios Ri{t) for the isospin eigenchannels / = and 2 are plotted in 
Fig. 10. To extract the energy shift SEj for each channel, we again employ a 
linear form Zj(l — SEjt) with the fitting range chosen to be 4 < i < 9. The 
values of the scattering lengths are listed in Table 5, together with those predicted 
by current algebra (42-43) but evaluated with the measured values of and 
/;r in Table 3. We quoted only the statistical errors for the scattering lengths 
Oq. Simulation results are consistent with current algebra and PCAC within 1-2 
standard deviations. For the 1 — 2 channel a similar agreement was previously 
found in Ref. [19]. 
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4.3 Question of quenched chiral divergences 

All of our results are obtained in quenched QCD. A problem with quenched 
calculations is the uncertainty in O(i^) term of the ratio Ri{t) due to lack of 
dynamical quark loops. We have bypassed this problem by looking only at the 
region of t where the ratio Ri{t) shows a hnear behavior. Also our data do not 
extend to the region of large t, where a curvature due to 0{t^) terms can be 
established, because of increase of errors. 

Another potential problem is that the tt-tt four-point function in quenched QCD 
could be affected by infrared divergences appearing in the chiral limit due to rj' 
loops[26]. These divergences originate from the double pole ml/{p'^ + ni^Y in the 
disconnected two quark loop amplitude of the rj' propagator, where ml — m^, —m^ 
represents the mass-squared splitting between the rj' and pseudoscalar octet mesons 
in the full theory. For the tt-tt scattering amplitude the one-loop diagram in chiral 
perturbation theory formed by two double-pole propagators yields a divergent 
imaginary part at threshold in the s channel [27]. This implies that the vacuum 
amplitude By it) would be ill-behaved as t increases. A rapid loss of signal seen in 
Fig. 4 and Fig. 9 may be partly related to this point. 

The same diagram viewed in the t channel gives a contribution of the form. 



to the s-wave scattering length, which diverges as — > 0. In terms of the energy 
shift in the direct amplitude R^{t), this contribution translates to 



with L the spatial size. 

Recently we have calculated [13] the parameter mo for quenched QCD with 
Wilson quark action using the same set of configurations on a 12^ x 20 lattice as 
are employed in the present work. At the hopping parameter K — 0.164 where 




(46) 




(47) 
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our scattering length calculation is made, we obtained mo/ ^jNf = 0.147(8). Com- 
bining this result with those of m.,^ and /jr in Table 3 we find SEd = 2.7 x 10~^ 
for the expected magnitude of energy shift for the direct amplitude for the Wilson 
case, which is far too small to be detectable in our data for the ratio (t) . 

The singular contribution (47) diverges in the chiral limit. In order to see at 
what quark mass the contribution seriously begins to affect the chiral behavior we 
recall[13] that the values for mo we obtained can be fitted very well with a linear 
function of the quark mass = {1/K — l/Kc)/2; 

0.299(14) - 1.57(19)mg. (48) 



/Nf 

For m^, mp and /^r we find from the values reported in Ref. [24] that 

ml = 2.72(10)mg, (49) 

mp = 0.5361(92) + 1.61(14)mq, (50) 

U - 0.0680(28) + 0.338(43)mg. (51) 

Combining (46-51) one can estimate the fractional change of the / scattering 
length relative to the current algebra value ao = 7m7r/(327r/^) expected for the 
Wilson action at P — 5.7. The result is plotted by a sohd line in Fig. 11 as 
a function of m.,^/mp. We see that the contribution of the singular term would 
become significant only for realistically small values of quark masses, which are 
beyond our computing resources. 

A similar estimate for the Kogut-Susskind case, strictly speaking, requires re- 
sults for mo for that action which is not available. One may, however, take over 
(48) as a guide. Employing the data of Ref. [28] for the Kogut-Susskind quark 
action at (3 — 5.7, which yield, 

ml = 0.004(1) + 7.96(5)mg, (52) 
mp = 0.782(86) + 10.0(6.7)m5, (53) 
= 0.118(7) + 0.97m5, (54) 
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one finds the dotted line in Fig. 11 for the fractional change 5ao/ao for the Kogut- 
Susskind case. Again the singular contribution is very small Sao/ao ~ 0.001 at the 
point nig — 0.01 {111^^/ nip — 0.33) of the simulation. For the energy shift we find 
5E = 3.8 X 10~^, which may be compared with the value 5E = 1.2(4.0) x 10~^ 
extracted from the slope of the direct amplitude R^{t). 

The estimate above indicates that the spurious infrared divergence and possible 
failure of the quenched approximation do not become manifest in tt-tt four-point 
functions unless the simulation is made with a quark mass much smaller than the 
value being taken in current quenched QCD studies. Of course, the problem of 
unphysical singularity does not arise if calculations are made on full QCD gauge 
configurations. 

4.4 Summary of tt-tt scattering lengths 

Our results for tt-tt scattering lengths are summarized in Fig. 12 in terms of 
the dimensionless ratio 327r/^ao/m7r together with those of Ref. [18, 19] for 1 = 2. 
We observe an agreement of lattice results with current algebra predictions up to 
quite heavy quark masses {ni-,^/nip ^ 0.7 — 0.8) for both 7 = and 2 channels. In 
detail, however, the Wilson result for the / = scattering length obtained at a 
heavy quark mass is somewhat smaller than the current algebra value, while the 
Kogut-Susskind results at a small quark mass are larger. It would be an interesting 
problem to repeat the Kogut-Susskind simulation for larger quark masses in order 
to see if the trend seen for the Wilson action is reproduced. For the 7 = 2 channel 
we do not observe such a difference between the two actions. 

Aside from the complications for the Kogut-Susskind case due to breaking of 
flavor symmetry discussed in Sec. 4.1.1, possible sources of systematic errors in our 
results are the uncertainties in higher order 1/L terms of (1) and scaling violations 
due to a fairly large lattice spacing of our simulation (a^^ ^ 1 — 1.5 GeV at (3 = 5.7 
determined from the p meson mass [2 3, 24]). For the 7 = 2 channel the previous 
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work[18, 19] checked these points. In particular, a comparison of results at P — 5.7 
and 6.0 indicates that finite lattice spacing errors are not large, at least for the 
Kogut-Susskind case, for the 1 — 2 channel (compare open squares {P — 5.7) and 
open diamonds (P ~ 6.0) in Fig. 12). Whether this applies also to the 7 = 
channel should be examined in future work. 

5 TT-N, K-N and K-N scattering lengths 

5.1 t:-N scattering lengths 

We calculate the t^-N scattering lengths from t^-N four-point functions. Since 
the wall source technique without gauge fixing does not yield good signals for the 
nucleon, we use the Coulomb gauge fixing at the t — Q time slice for the nucleon 
source. We employ the Wilson quark action at hopping parameters K — 0.164 
and K = 0.1665, using 60 and 30 gauge configurations, respectively. 

In Fig. 13 we show the ratio Ri{t) of tt-N amplitude for the isospin eigenchan- 
nels 1 — 1/2 and 3/2 at the two values of the hopping parameter. The data are 
fitted with the linear form Ri{t) = Zj{l — 6Ejt) over 4 < t < 9, and the results for 
the scattering length are listed in Table 6. They are compared with predictions of 
current algebra and PCAC (values in brackets in Table 6), which are given by 



where /It^n — mT^niN /{ttIt^ + 'mN) is the reduced mass and the right-hand sides are 
estimated with the measured values of 171^^,171 n and /tt in Table 3. The results at 
K = 0.164 are consistent with the current algebra prediction, in spite of the fact 
that simulations are made with quite heavy quarks corresponding to m^/m^ = 
0.74. The quality of data deteriorates for a fighter quark oi K — 0.1665; more 
statistics are clearly needed as Fig. 13 indicates. 



ao(/ = l/2) = 4 



ao{I = 3/2) 




(55) 



(56) 
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Let us define (t) to be the contribution of X type diagrams {X — D,C,R 
and CR; see Fig. 2) to crossing even (+) and odd (— ) cliannels, summed over 
fermion contractions having the quark hne topology of the X diagram. In terms 
of the isospin eigen amphtudes Ri{t), the crossing even and odd amphtudes are 
given by 

R+{t) = ^Ri=i/2{t) + lRi=3/2{t), (57) 
R4t) = lRi=i/2{t)-lRi=3/2{t). (58) 

The ratios R^{t) are plotted in Fig. 14 as a function of t ior K — 0.164. We recall 
that the crossed (C) and rectangular (R) amplitudes are governed by current 
algebra; the soft pion theorem predicts that the crossing even amplitudes R^{t) 
and R'^{t) should show a reversal of sign in the slope with respect to t. The direct 
(D) and crossed rectangular (CR) amplitudes do not appear in current algebra to 
leading order in m.,^, and therefore we expect that R±{t) and i?^^(t) are fiat. Our 
results in Fig. 14 are roughly consistent with these expectations except that the t 
dependence is not very pronounced for the rectangular diagram compared to the 
crossed diagram and that the direct amplitude exhibits a positive curvature. Both 
these features are similar to what appeared in the tt-tt case with the Wilson quark 
action, perhaps originating from heavy quark mass employed {mT^/mp — 0.74). 

5.2 K-N and K-N scattering lengths 

The prediction of current algebra and PCAC for SU(3) symmetry is given by 



4\I = Q) = 0, (59) 

ar(^ = l) - (60) 

ar(/ = 0) = 1^^, (61) 

ar(/ = l) = (62) 
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with fiKN — f^KiT^N I {triK + i^n) being the reduced mass. Looking at Table 1 
we observe that the results for the KN channel are not unreasonable compared 
with experiment, although the I — 1 scattering length is off by a factor of two 
probably due to a long extrapolation to the soft K meson limit. On the other hand, 
comparison of the predictions (61-62) for the K-N channel with the experiment 
requires a caution: ttA and ttE channels are already open at the KN threshold. It 
is interesting to note that the prediction (62) for the I —1 KN scattering length 
still agrees with the real part of the experimental value. The fact that the real 
part of the 7 = KN scattering length predicted by (61) has a sign opposite to 
that of the experiment is clearly due to the presence of A(1405). We point out that 
the lattice results reported below can be compared better with the current algebra 
values, since the present lattice calculation with the quenched approximation is 
basically a single channel calculation without taking into account an opening of 
the channels below the threshold. 

Our results for K-N and K-N scattering lengths are obtained at K — 0.164, 
employing the same hopping parameter for the strange and up-down quark masses. 
The contribution of the direct and crossed diagrams to the K-N amplitude in the 
7 = and 1 isospin eigenchannels is plotted in Fig. 15. For the / = channel 
(filled circles) a small positive slope of the direct amphtude is almost canceled out 
by the negative slope of the crossed one. The direct and rectangular contribution 
to the K-N amplitude for both isospin eigenchannels are presented in Fig. 16. 

The ratio Ri{t) for the isospin eigenchannels is shown for the K-N and K-N 
cases in Fig. 17 and 18 respectively. We made a linear fit Riit) = Zj{l — 6Ejt) 
over 4 <t <9 with the results given in Table 7, together with the current algebra 
values evaluated with the measured value of rriT^^mj^ and /^r in Table 3. We see 
that all signs agree with those of the current algebra results. For the 1 = KN 
scattering length for which the current algebra value is vanishing, the lattice result 
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also yields a smaller value. For other cases, the agreement is up to a factor of two. 

6 N-N scattering lengths 

Calculation of N-N scattering lengths poses a number of challenging obstacles 
not encountered for the case of tt-tt and tt-N scattering. First of all the experi- 
mental scattering lengths are quite large, being of order lOfm (see Table 1). This 
means that lattice sizes much larger than 2 — 3fm, which are accessible in current 
numerical simulations, will be needed to suppress 0{L~^) corrections neglected in 
the Liischer's formula (1). Secondly extraction of the scattering length in the spin 
triplet channel requires a calculation of the lowest scattering state orthogonal to 
the ground state, which is the deuteron bound state. In fact the negative sign 
of the s-wave spin triplet scattering length is a consequence of the existence of 
the deuteron as follows from Levinson's theorem. Finally statistical fluctuations 
in N-N four-point functions are expected to grow more rapidly than in the tt-tt 
and TT-A^ cases toward large time separations and small quark masses as we shall 
discuss in detail in Sec. 7 below. 

A possible strategy in this situation is to start from the region of heavy quark 
masses, where simulations are easier, and subsequently reduce the quark mass. 
To follow this approach, we first examine the behavior of nucleon-nucleon forces 
toward larger quark masses through a phenomenological model of one-boson ex- 
change potentials. We then present our lattice results carried out with heavy 
quarks corresponding to rriT^/mp ^ 0.74 — 0.95. 

6.1 Phenomenological considerations 

Low energy N-N scattering data are usually described by phenomenological 
models of one-boson exchange potentials [3]. At large nucleon-nucleon separations 
r ^ 2fm, the potential is dominated by one-pion exchange. At intermediate dis- 
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tances 2 ^ r ^ Ifm, the two-pion exchange provides the dominant attraction, 
which is modelled by an isoscalar scalar particle a of mass ~ 500MeV. Heavier 
bosons such as p and lo also become important in this range of distance. The 
potential turns into a hard core at r ~ 0.5fm whose dynamical details are not 
understood well. 

The formation of bound states in the triplet ^5"! channel that emerges from 
such a potential is subtle. The central potential, which has a depth of about 
— 50MeV at r ~ Ifm both in the triplet and singlet ^5*0 channels, is insufficient 
for bound state formation. For the triplet channel, however, the tensor potential 
provides an additional attraction of a similar depth. This leads to the binding of 
the deuteron with a small binding energy of 2.22457MeV with a sizable mixing of 
D wave. 

Let us consider what happens if the quark mass is increased from the physical 
value. The mesons most affected by the increase is the pion and a representing 
the two-pion exchange, whose masses increase as ^/m^. Thus the attractive po- 
tential provided by these states becomes rapidly shorter ranged (while the hard 
core mostly given by lo changes little), and the two nucleons have to come closer 
to remain in the potential well. One can imagine that the resulting increase of 
kinetic energy may easily overwhelm the small deuteron binding energy, leading 
to unbinding of the deuteron state. 

We have examined whether this actually takes place by employing the one- 
boson exchange model of Ref. [29, 30] and varying the values of ttItt and m„ ac- 
cording to our lattice result for the slope rn^/mq. We find a divergence of the 
scattering length for the triplet channel taking place at rUq — 6.3MeV that signals 
unbinding of the deuteron. The fractional increase of quark mass from the value 
4.9MeV corresponding to the physical point is only 30%. We have repeated the 
analysis in several alternative ways, e.g., varying other meson masses and/or the 
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nucleon mass assuming the relation mn — An + BHrng, and also including a varia- 
tion of the meson-nucleon coupling constants. We have found that the qualitative 
features are not modified: the deuteron most probably ceases to exist as a bound 
state when the quark mass is increased by more than 30-40% from the physical 
value. Another common feature is that the scattering lengths decrease to a value 
of the order of Ifm, when quark mass is increased so that rriT^/mp ~ 0.7 — 0.9. It is 
reasonable to expect that the scattering length takes a value similar to the hadron 
size in the absence of bound-state effects. 

6.2 Lattice results 

Our lattice study is carried out with the Wilson quark action. In the calculation 
of the nucleon four-point function gauge configurations are fixed to the Coulomb 
gauge over all space-time to enhance signals. Projecting out spin singlet and triplet 
combinations of the two nucleon system is made by non-relativistically combining 
the upper Dirac components of the nucleon operator as explained in Sec. 2.2. Other 
technical features are similar to the tt-tt and t:-N cases. 

In Fig. 19 we show the individual contributions of the direct and crossed am- 
plitudes to the ratio R{t) for the triplet ^5*1 (filled circles) and singlet ^5*0 (open 
triangles) channels obtained at = 0.160. We find that the direct amplitudes for 
the two channels are virtually the same. On the other hand, the crossed amplitude 
in the spin triplet channel increases in contrast to a fiat behavior in the spin singlet 
channel. 

In Fig. 20 we show R{t) for the spin singlet and triplet channels aX K — 0.160, 
which corresponds to 171^^ /nip = 0.85. A clear signal with a positive slope is 
observed for both channels, which means attraction {SEnn — Enn — ^uin < 0). 
Similar results are obtained at two other values of the hopping parameter K — 
0.150(m^/mp = 0.95) and 0.164(7727^/^^ = 0.74). We extract the energy shift 
5Enn by fitting R{t) to a hnear form R{t) — Z{1 — SEnnI). The fitting range 
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is chosen to be 4 < i < 9 for X 0.150 and 0.160. For the case oi K ^ 0.164, 
however, we used 2 < t < 6 due to poor quality of our data. The fitted values of 
SEj^N are quite small (~ 0.01), justifying the use of a linear function instead of an 
exponential. 

From phenomenological considerations in Sec. 6.1 we expect that the deuteron 
is not a bound state at a heavy quark mass where our simulations are made. We 
then extract the scattering lengths through (1) for both spin singlet and triplet 
channels. The results in lattice units are tabulated in Table 8. We should remark 
that the scattering lengths we found are large enough to warrant a calculation 
with yet a larger lattice size, even though we used quite a large size of 20^: the 
three terms in (1) are comparable in magnitude. A finite- size analysis with larger 
lattice sizes is clearly necessary. 

In Fig. 21 we compare the results for the N-N scattering lengths with those for 
Ti-n and vr-A^ obtained with the Wilson quark action. Conversion to physical units 
is made using a — 0. 137(2) fm determined from the p meson mass. It is apparent 
that the N-N scattering lengths are substantially larger than the n-N and tt-tt 
scattering lengths aheady for a heavy quark corresponding to mT^/mp ~ 0.74. 
Also noteworthy is the trend, albeit with sizable errors, that the values for the 
spin triplet channel are larger than those for the singlet channel, indicating a 
stronger attraction in the triplet channel. This is consistent with the existence of 
the deuteron bound state for physical quark mass. 

In physical units our results correspond to ao{NN) 1.0 — 1.5fm. These val- 
ues are small compared to the experimental value of order lOfm. However, the 
large experimental scattering lengths are a reflection of the fact that the s-wave 
nucleon-nucleon system is either marginally bound (triplet channel) or very close to 
a bound-state formation (singlet channel), which would not be the case for heavy 
quarks studied in our simulation. Therefore, our results are not unreasonable. We 
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expect that lattice results will exhibit an increase as the quark mass is decreased. 
For the triplet channel, in particular, the scattering length should diverge at the 
point where the deuteron bound state is formed. According to the phenomenolog- 
ical consideration in the previous section, however, this will take place quite close 
to physical quark mass, which is not accessible in current lattice QCD simulations. 

7 Problem of statistical errors 

Our results for R{t) for various channels share the feature that statistical errors 
grow rapidly for large times, examples of which are seen in Fig. 5 and Fig. 10 for 
TT-TT, Fig. 13 for TT-N and Fig. 20 for N-N cases. This behavior can be understood 
from (41), which implies that the error SR^{t) of the ratio R^{t) for a diagram 
X increases as SR^{t) oc exp(Q;x^) with the exponent ax given in Table 9. The t 
dependence of the measured errors agree very well with this formula. Some typical 
examples are shown in Fig. 22. 

This consideration clarifies why the vacuum amplitude for the tt-tt scattering 
rapidly loses signal as t increases: the error has the largest exponent ay = 2m7r 
among the four amplitudes contributing to the tt-tt four-point function. We also 
expect that the rate of growth of errors decreases for the tt-tt case toward small 
quark mass, since the exponent is governed by the pion mass. Exploring the light 
quark mass region is more difficult for tt-A'^ and N-N scattering: the exponent for 
these cases is a difference of the nucleon mass and pion mass with some coefficient, 
which depends on the channel, and the exponent stays non-vanishing or even 
increases toward the chiral limit. We illustrate this point in Fig. 23 where we plot 
the largest exponent, 2-171^^ for tt-tt, mjv — ?77.^/2 for tt-N and 2m^ — 3m^ for N-N 
function of m.„/mp where we used the fits (49-50) for and nip and 

mjv = 0.741(21) + 3.80(31)mg (63) 

for niN obtained from the hadron mass data in Ref. [24]. We observe that the 
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problem will be particularly severe for the A'^-A'^ four-point function for which the 
exponent begins to increase for mT^/mp ^ 0.5. 

8 Conclusions 

In this work we have calculated hadron scattering lengths in a variety of chan- 
nels at (5 — 5.7 in the quenched approximation. We have demonstrated that 
Liischer's formula that correlates scattering lengths with the energy of the two- 
hadron system confined in a finite box can be used to calculate any class of di- 
agrams of scattering amplitudes, when implemented with a modified wall source 
method proposed by the present authors. The results arc physically sensible and 
encouraging, being in a reasonable agreement with current algebra values and the 
experiment for tt-tt and t^-N cases; for N-N scattering the resulting scattering 
lengths are significantly larger compared to the first two cases even for heavy 
quarks. 

In the present study dynamical effects of sea quarks are not taken into account. 
However, it is straightforward to apply our method once configurations are gener- 
ated with dynamical quarks. Worrisome toward a more realistic calculation is the 
exponential growth of errors for large times discussed in Sec. 7. The problem is 
particularly acute for N-N scattering, for which the rate is numerically large and 
is expected to increase as the chiral limit is approached. The anticipated large 
N-N scattering lengths, which henceforth requires a large lattice size, makes the 
problem further difficult. Advance in our understanding of the low energy nucleon- 
nucleon scattering from the first principle has to await not only future progress of 
computing power but also further innovation of calculational techniques. 
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Appendix 

In this Appendix we enumerate weights of various quark contractions con- 
tributing to the TT-N, K-N and ^-A^ four-point functions. Let us consider the 
meson-baryon operators given by 

0{x,y) = eojbc{^qiC-^lf,qt)(ll{x)xq^lf,qi{y). (64) 
0\x\y') = ea'b'c4{i{C-'75V'qt'){x')x-q,^,q-,{y'). (65) 

where subscripts with tilde are assigned to quark fields q and those without tilde to 
anti-quark fields q. Quark contractions for the operator product 0{x,y)0\x' ,y') 
may be specified by permutations from the set {1, 2, 3, 4, 5} for quark fields to the 
set {1,2,3,4,5} for anti-quark fields. For example, the contraction with crossed 
rectangular topology shown in Fig. 24 corresponds to 

12 3 4 5 
4 5 2 3 1. 

There are 6, 18, 18 and 36 inequivalent contractions for the direct, crossed, rect- 
angular and crossed rectangular diagrams. In Table 10 - 13 we summarize the 
weights of the corresponding amplitudes for tt-N, K-N and K-N four-point func- 
tions projected to isospin eigenchannels. Sign factors arising from Fermi statistics 
are not included in the weight. For definitions of meson and nucleon operators, 
see Sec. 2.2. 
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Tables 



Table 1: Experimental s-wave scattering lengths[l] and current algebra predictions 
in units of fm. 









experiment 


current algebra 


TT-TT 


I = 





+0.37(7) 


+0.222 




I = 


2 


-0.040(17) 


-0.0635 


TT-N 


I = 


1/2 


+0.245(4) 


+0.221 




7 = 


3/2 


-0.143(6) 


-0.111 


N-N 






-5.432(5) 










+20.1(4) 




K-N 


I = 





-0.0075 







I = 


1 


-0.225 


-0.399 


K-N 


1 = 





-1.16 + i0.49 


+0.598 




I = 


1 


+0.17 + /0. 11 


+0.199 



Table 2: Parameters of simulation. All runs are made at l3 = 5.7 in quenched 
QCD. 





action 


quark mass 


lattice size 


# conf. 




gauge fixing 


TT-TT 


KS 


ruq = 0.01 


8^ X 20 


400 




none 




KS 


rUg = 0.01 


12^ X 20 


160 




none 




KS 


nig = 0.01 


12^ X 20 


60 




Coulomb 




Wilson 


K = 0.164 


12''^ X 20 


70 




none 


TT-N 


Wilson 


K = 0.164 


12''^ X 20 


60 


t 


= Coulomb 




Wilson 


K = 0.1665 


12^ X 20 


30 


t 


= Coulomb 


K-N 


Wilson 


K = 0.164 


12^ X 20 


60 


t 


= Coulomb 


K-N 


Wilson 


K = 0.164 


12=^ X 20 


60 


t 


= Coulomb 


N-N 


Wilson 


K = 0.150 


20^ X 20 


20 




Coulomb 




Wilson 


K = 0.160 


20''^ X 20 


30 




Coulomb 




Wilson 


K = 0.164 


20^ X 20 


20 




Coulomb 



Table 3: Hadron masses and pion decay constant at /9 = 5.7 in quenched QCD. 
Also listed for comparison are results on a larger lattice from Refs. [23, 24]. For 
the Wilson result of Ref . [24] a hnear interpolation in 1/K \s made when necessary 
(marked by *). 



action quark mass 


lattice size 




rriN 




KS ruq = 


-- 0.01 


12^ X 20 


0.290(3) 




0.132(3) 






24'^ X 32[23] 0.2876(7) 0.883(58) 


1.454(26) 




Wilson K = 


0.150 


12^ X 20 


1.0758(51) 1.1302(66) 




0.1900(12) 






20=^ X 20 




1.788(11) 




K = 


0.160 


12^ X 20 


0.6876(31) 0.8053(39) 


1.2957(79) 


0.1268(14) 






20=^ X 20 




1.302(13) 








24=^ X 32 [24] 


0.6887(11) 0.8062(18) 


1.3068(51) 


0.1249(10) 




0.164 


12^ X 20 


0.5080(37) 0.6865(53) 


1.080(10) 


0.1009(7) 






20^ X 20 




1.093(20) 








24^ X 32 [24] 


0.5027(13)* 0.6886(28)* 


1.0869(64)* 


0.1000(10)* 




0.1665 12^ X 20 


0.3663(44) 0.6086(75) 


0.926(13) 


0.0832(14) 






243 X 32 [24] 


0.3656(17)* 0.6165(45)* 0.927(12)* 


0.0848(16)* 



Table 4: tt-tt scattering lengths in lattice units for the Kogut-Susskind quark 
action at rUg = 0.01. Numbers in square brackets are current algebra predictions 
evaluated with the measured values of TOtt and listed in Table 3. 



gauge 7 = 1 — 2 

size fixing oq 5Ei Zj oq 5Ei Zj 

8^ X 20 iw^e 4.89(53) -0.69(21) 0.160(40) -0.374(14) 0.0367(16) 0.897(12) 
12^ X 20 none 1.57(25) -0.0291(37) 0.807(15) -0.301(28) 0.00813(82) 0.955(5) 
Coulomb 1.73(27) -0.0316(40) 0.800(17) -0.326(35) 0.0089(10) 0.948(7) 
[1.16(5)] [-0.331(15)] 
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Tabic 5: tt-tt scattering lengths in lattice units for the Wilson quark action at the 
hopping parameter K = 0.164 on a 12^ x 20 lattice. Numbers in square brackets 
are current algebra predictions evaluated with the measured values of and f^^ 
listed in Table 3. 





7 = 


7 = 2 




ao 


SEi Zi 


ao SEi 


Zi 


3.02(17) 


-0.0297(19) 0.903(7) 


-0.924(40) 0.0166(9) 


1.027(5) 


[3.47(5)] 




[-0.993(16)] 





Table 6: ir-N scattering lengths in lattice units for the Wilson quark action on 
a 12^ X 20 lattice. Numbers in square brackets are current algebra predictions 
evaluated with the measured values of TOtt, mjv and f^^ listed in Table 3. 







7 = 1/2 




7 = 3/2 


K 


ao 


SEi Zi 


ao 


SEi Zi 


0.164 


3.04(66) 


-0.0219(54) 0.916(22) 


-1.10(20) 


0.0151(35) 0.975(18) 




[2.701(41)] 




[-1.350(20)] 




0.1665 


-0.70(58) 


0.011(11) 1.021(52) 


-1.31(22) 


0.0243(55) 0.964(31) 




[3.02(11)] 




[-1.509(53)] 





Table 7: K-N and K-N scattering lengths in lattice units for the Wilson quark 
action at the hopping parameter K = 0.164 on a 12'^ x 20 lattice. Numbers in 
square brackets are current algebra predictions evaluated with the measured values 
of mTT, mjv and /^r listed in Table 3. 





7 = 




7 = 1 


ao 


SEi Zi 


ao 


SEi Zi 


K-N 0.55(47) 


-0.0051(38) 0.987(19) 


-1.56(13) 


0.0240(27) 1.012(16) 


[0] 




[-2.701(41)] 




K-N 4.64(37) 


-0.0415(64) 0.881(23) 


2.63(64) 


-0.0188(45) 0.941(20) 


[4.051(61)] 




[1.350(20)] 
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Table 8: N-N scattering lengths in lattice units for the Wilson quark action on a 
20^ lattice. 



^1 

K ao 5^ Z oo 5^ Z 

0.150 10.8(1.2) -0.0126(39) 1.012(19) 9.2(1.3) -0.0085(27) 0.987(13) 
0.160 9.0(1.6) -0.0109(45) 1.048(25) 7.3(1.9) -0.0072(33) 1.007(18) 
0.164 10.8(9) -0.0207(48) 1.021(12) 8.0(1.1) -0.0102(30) 0.997(7) 



Table 9: Exponent a governing the growth of error of the ratio R^{t) for a daigram 
of type X. 



X D C R V CR 

TT-TT 2m7r 

TT-N rriN — |?77.7r ttin — Itott uin — \m-^ rriN — ^ttIt^ 

N-N 2mjv — 3m^ 2mjv — 3m^ 



Table 10: Weights of direct (D) diagrams for tt-N, K-N and K-N four-point func- 
tions projected to isospin eigenchannels. 







TT- 


■N 




K-N 




K-N 




diagram 


contraction 


7 = 1/2 


I 


= 3/2 


7 = 7 


= 1 


7 = 7 


= 1 


Dl 


12345 


2 




2 


2 


2 


2 


2 


D2 


13245 


1 




1 


1 


1 


1 


1 


D3 


21345 


-2 




-2 


-2 


-2 


-2 


-2 


D4 


31245 


-1 




-1 


-1 


-1 


-1 


-1 


D5 


32145 


1 




1 


1 


1 


1 


1 


D6 


23145 


-1 




-1 


-1 


-1 


-1 


-1 
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Table 11: Weights of crossed(C) diagrams for n-N, K-N and K-N four-point 
functions projected to isospin eigenchannels. 



diagram 


contraction 


TT-N 

1=1/2 I 


= 3/2 


K-N 
1 = I 


= 1 


K-N 
1 = I 


= 1 


CI 


42315 


1 


1 


1 


1 








C2 


43215 


3/2 





2 











C3 


24315 


-1 


-1 


-1 


-1 








C4 


34215 


-3/2 





-2 











C5 


23415 


1/2 


-1 


1 


-1 








C6 


32415 


-1/2 


1 


-1 


1 








C7 


41325 


-1 


-1 


-1 


-1 








C8 


43125 


-3/2 





-2 











C9 


14325 


1 


1 


1 


1 








CIO 


34125 


3/2 





2 











Cll 


13425 


-1/2 


1 


-1 


1 








C12 


31425 


1/2 


-1 


1 


-1 








C13 


41235 


1/2 


-1 


1 


-1 








C14 


42135 


-1/2 


1 


-1 


1 








C15 


14235 


-1/2 


1 


-1 


1 








C16 


24135 


1/2 


-1 


1 


-1 








C17 


12435 


-1 


2 


-2 


2 








C18 


21435 


1 


-2 


2 


-2 
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Table 12: Weights of rectangular (R) diagrams for t^-N, K-N and K-N four-point 
functions projected to isospin eigenchannels. 



diagram 


contraction 


TT-N 

1=1/2 I 


= 3/2 


K-N 
1 = I 


= 1 


K-N 
1 = I 


= 1 


Rl 


52341 


1 


1 








1 


1 


R2 


53241 


-1/2 


1 








-1 


1 


R3 


25341 


-1 


-1 








-1 


-1 


R4 


35241 


1/2 


-1 








1 


-1 


R5 


23541 


-3/2 











-2 





R6 


32541 


3/2 











2 





R7 


51342 


-1 


-1 








-1 


-1 


R8 


53142 


1/2 


-1 








1 


-1 


R9 


15342 


1 


1 








1 


1 


RIO 


35142 


-1/2 


1 








-1 


1 


Rll 


13542 


3/2 











2 





R12 


31542 


-3/2 











-2 





R13 


51243 


-3/2 











-2 





R14 


52143 


3/2 











2 





R15 


15243 


3/2 











2 





R16 


25143 


-3/2 











-2 





R17 


12543 


3 











4 





R18 


21543 


-3 











-4 
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Table 13: Weights of crossed rectangular(CR) diagrams for n-N, K-N and K-N 
four-point functions projected to isospin eigenchannels. 







TT-iV 




K-N 




K-N 




diagram 


contraction 


I — 1/2 I 


— 3/2 


I — I 


= 1 


I — I 


= 1 


CRl 


52431 


-1/2 


1 














CR2 


54231 


-1/2 


1 














CR3 


53421 


-1/2 


1 














CR4 


54321 


1 


1 














CR5 


25431 


1/2 


-1 














CR6 


45231 


1/2 


-1 














CR7 


35421 


1/2 


-1 














CR8 


45321 


-1 


-1 














CR9 


24531 




















CRIO 


42531 




















CRll 


34521 


3/2 

















CR12 


43521 


-3/2 

















CR13 


51432 


1/2 


—1 














CR14 


54132 


1/2 


-1 














CR15 


53412 


1/2 


-1 














CR16 


54312 


-1 


-1 














CR17 


15432 


-1/2 


1 














CR18 


45132 


-1/2 


1 














CR19 


35412 


-1/2 


1 














CR20 


45312 


1 


1 














CR21 


14532 




















CR22 


41532 




















CR23 


34512 


-3/2 

















CR24 


43512 


3/2 



















0i4zo 


U 


A 

u 


u 


A 

u 


A 

u 


A 

u 


CR26 


54123 


3/2 

















CR27 


52413 




















CR28 


54213 


-3/2 

















CR29 


15423 




















CR30 


45123 


-3/2 

















CR31 


25413 




















CR32 


45213 


3/2 

















CR33 


14523 


3/2 

















CR34 


41523 


-3/2 

















CR35 


24513 


-3/2 

















CR36 


42513 


3/2 
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Figure Captions 



Figure 1: Diagrams contributing to tt-tt four-point functions. Short bars represent 
wall sources. Open circles are sinks for local pion operator. 

Figure 2: Diagrams contributing to tt-N four-point functions. Short bars represent 
wall sources. Open circles are sinks for local pion or nucleon operator. Solid circle 
for the CR diagram means use of source method for calculating the product of 
quark propagators connected by it. 

Figure 3: Diagrams contributing to A^-A^ four-point functions. Short bars represent 
wall sources. Open circles are sinks for local nucleon operator. 

Figure 4: Individual ratios R^{t) {X = D,C,R and V) for tt-tt four-point function 
calculated with the method of wall source without gauge fixing for Kogut-Susskind 
quark action on a 12^ x 20 lattice at = 5.7 and rriq = 0.01. 

Figure 5: Ri{t) {1 = and 2) for tt-tt four-point function calculated without gauge 
fixing for the Kogut-Susskind quark action on a 12^ X 20 lattice at /3 = 5.7 and 
iTiq = 0.01. Solid lines arc linear fits for 4 < t < 9. 

Figure 6: Ri(t) {1 = and 2) for tt-tt four-point function calculated without gauge 
fixing for the Kogut-Susskind quark action on an 8^ x 20 lattice at /3 = 5.7 and 
niq = 0.01. Solid lines are linear fits for 4 < t < 9. 

Figure 7: Comparison of Ri{t) {1 = and 2) for tt-tt four-point function calculated 
in Coulomb gauge (open symbols) and without gauge fixing (filled symbols) for the 
Kogut-Susskind quark action for the same set of 60 configurations on a 12^ x 20 
lattice at /3 = 5.7 and ruq = 0.01. 

Figure 8: Comparison of single-elimination jackknife errors of Rj{t) {1 = and 
2) for TT-TT four-point function calculated in Coulomb gauge (open symbols) and 
without gauge fixing (filled symbols) for the Kogut-Susskind quark action for the 
same set of 60 configurations on a 12^ x 20 lattice at /3 = 5.7 and niq = 0.01. 

Figure 9: Individual ratios R^{t) {X = D,C,R and V) for tt-tt four-point function 
calclulated without gauge fixing for Wilson quark action on a 12^ x 20 lattice at 
13 = 5.7 and K = 0.164. 

Figure 10: Ri{t) {1 = and 2) for tt-tt four-point function calculated without 
gauge fixing for the Wilson quark action on a 12^ x 20 lattice at /3 = 5.7 and 
X = 0.164. Solid lines are linear fits for 4 < i < 9. 

Figure 11: Fractional change of the 7 = tt-tt scattering length due to in- 
frared divergences in quenched QCD normalized by the current algebra value 
flo = 7m^/ (32tt/^) at j3 = 5.7 plotted as a function of uiT^/mp. Solid line is for the 
Wilson action and dotted line for the Kogut-Susskind action. Arrows indicate the 
position where our calculation of tt-tt four-point functions is made. 
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Figure 12: 7 = and 2 s-wave n-n scattering lengths Oq. Filled and open symbols 
denote Wilson and Kogut-Susskind results. Triangles are for Coulomb gauge re- 
sults. Squares {(3 — 5.7) and diamonds (/3 = 6.0) for 7 = 2 are from Refs. [18, 19]. 
Dotted lines indicate predictions of current algebra. 

Figure 13: Ri{t) (7 = 1/2 and 3/2) for the n-N four-point function at (a) K = 
0.164 and (b) 0.1665. Solid hues are linear fits for 4 < i < 9. 

Figure 14: Individual ratios R±{t) {X = D,C,R and CR) for tt-N four-point 
function projected to crossing even (filled circles) and odd (open circles) channels. 

Figure 15: Individual ratios Rf{t) {X = D and C) of the K-N amplitude for 
isospin eigenchannels. Filled and open symbols denote 7 = and 1 results. 

Figure 16: Individual ratios Rf{t) {X = D and 7?) of the K-N amplitude for 
isospin eigenchannels. Filled and open symbols denote 7 = and 1 results. 

Figure 17: Ri{t) (7 = and 1) for the K-N four-point function on a 12^ x 20 
lattice at /3 = 5.7 and K = 0.164. Solid lines are linear fits for 4 < t < 9. 

Figure 18: Ri{t) (7 = and 1) for the ^-A^ four-point function on a 12^ x 20 
lattice at P — 5.7 and = 0.164. Solid hues are linear fits for 4 < i < 9. 

Figure 19: Individual contributions of direct and crossed amphtudes to the ratio 
R{t) for triplet ^y^i (filled symbols) and singlet ^5'o(open symbols) channels. 

Figure 20: R{t) for N-N four point function on a 20^ lattice at /3 = 5.7 and 
K = 0.160. Solid hues are linear fits for 4 < i < 9. 

Figure 21: N-N scattering lengths in units of fm as compared to n-N and n-ir 
scattering lenghts in quenched QCD at P — 5.7 calculated with the Wilson quark 
action. Conversion to physical units is made with a — 0.137(2)fm determined from 
the p meson mass. 

Figure 22: Error SR^{t) of ratio R^{t) as a function of i at /3 = 5.7 in qenched 
QCD. Use of Kogut-Susslind or Wilson quark action is indicated in parentheses. 
Solid lines indicate expected slope calculated with measured values of rriN and m^. 
(a) vr-TT at = 0.01, (b) vt-tt at 7s: = 0.164, (c) tt-N at K = 0.164, (d) N-N in 
the triplet '^^i channel at K = 0.160. 

Figure 23: Largest exponent a governing the errors of ratio SR{t) oc e"* for tt-tt 
(dotted line), tt-N (broken line) and N-N (solid line) cases. 

Figure 24: Example of quark contractions contributing to tt-N four-point function. 
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